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We present results of a search for the X(3872) in B'^ X{3S72)K'^,X{3872) J/i/'Tr+Tr", 
improved measurements of B" X{3872)K^ , and a study of the J/ipn'^n~ mass region above 
the X(3872). We use 232 milhon BB pairs collected at the r(4S) resonance with the BABAR 
detector at the PEP-II e"'"e~ asymmetric-energy storage rings. The results include the 90% con- 
fidence interval 1.34 x 10"^ < B{B° -* X{3872)K'',X -> J/t/^tt+tt") < 10.3 x 10"*^ and the 
branching fraction 6(5" ^ X{3872)K-,X J /■>Ijti+-r-) = (10.1 ± 2.5 ± 1.0) x 10"^ We 
observe a (2.7 ± 1.3 ± 0.2)MeV/c^ mass difference of the X(3872) produced in the two decay 
modes. Furthermore, we search for the y(4260) in B decays and set the 95% C.L. upper limit 
B{B- ^(4260)7^-, r (4260) ^ J/iPtv+tt ) < 2.9 x 10'^ . 

PACS numbers: 13.25.Hw, 14.40.Gx, 12.39.Mk 



The X(3872) was first observed in the charged 
S-meson decay Q -B" X{3872)K-, X{3872) 
J/ipn~^TT~ by the Belle Collaboration [al- It has been 
confirmed by the BABAR Collaboration |3 and observed 
inclusively in the same final state by the CDF and DO col- 
laborations This narrovi^-v^ridth particle has a mass 
very near the Z)"/)*" threshold and decays into final 
states containing charmonium (J/ip). The most plausi- 
ble interpretation was a l'^D2 or 1^D2 cc state which 
would be narrow since it would be forbidden to decay 
into open charm DD states. However, these candidates 
should have large radiative transitions into Xc states that 
have not been observed 0- Recent studies from Belle 
that combine angular and kinematic properties of the 
TT+TT" mass, strongly favor a J^^ = 1++ state 0- Other 
explanations include 2^ Pi cc (1^ ) or 2^ Pi 1+"*" states 
that should be narrow, but are predicted to be about 
lOOMeV/c^ higher than the X{3872) mass and are not 
expected to have a large decay rate into J/tJjt:?: final 
states 0. Hence, the X{3872) appears not to be a simple 
quark model qq meson state. 

Many explanations have been proposed for the na- 
ture of the X(3872). Recent interpretations include the 
diquark-antidiquark model and the S-wave D^D*^ 
molecule model . The diquark-antidiquark model pre- 
dicts a spectrum of J = 0, 1, 2 particles and identifies the 
X(3872) as its l"*""*" member state with the two quark 
combinations Xu = [cu] [cu] and Xd = [cd\ [c3] with a 
mass difference m{Xd) ~ m{Xu) w (7 ± 2)M(N/c^. In 
addition, these two states could form mixed states that 
are produced in both charged and neutral _B-meson de- 
cays with different masses and rates depending on the 
mixing angle. A search for the predicted charged partner 
of the X(3872) has been addressed in a previous analy- 
sis with an upper limit that is still consistent with the 
model [13. The D^D*° molecule model interprets the 
X(3872) as a loosely bound D°D*° S-wave state t_hat is 
produced in weak decays of the i3-meson into D^D*^K. 



In this picture, the S-wave molecule must form a = 1"*" 
state. From factorization, heavy-quark symmetry, and 
isospin symmetry, the decay B^ X(3872)i^° is pre- 
dicted to be suppressed by an order of magnitude relative 
to B- X{3872)K- To investigate these predic- 
tions, we present in this letter a study of the neutral mode 
^ X{3872)K°, X{3872) J/i/^Tr+Tr" and we analyze 
B^ J/ipTT^TT^ decays with increased statistics to 
obtain improved measurements of X(3872) J/ipir'^TT^ . 
In addition, we examine the higher J/'07r"'"7r~ invariant 
mass region to search for a structure recently observed 
in initial state radiation (ISR) events ^3 • 

The data were collected with the BABAR detector at 
the PEP-II asymmetric-energy e~^e^ storage rings on the 
T{4:S) resonance. The integrated luminosity of the data 
used in this analysis is 211 fb^^; this corresponds to the 
production of (232 ± 3) x 10^ BB pairs. 

The BABAR detector is described in detail else- 
where 0. Charged-particle trajectories are measured 
by a combination of a five-layer silicon vertex tracker 
(SVT) and a 40-layer drift chamber (DCH) in a 1.5-T 
solenoidal magnetic field. For charged-particle identi- 
fication, we combine information from a ring-imaging 
Chcrcnkov detector (DIRC) and energy-loss measure- 
ments provided by the SVT and the DCH. Photons 
and electrons are detected in a CsI(Tl) electromagnetic 
calorimeter (EMC). Penetrating muons are identified by 
resistive-plate chambers in the instrumented magnetic 
fiux return (IFR). 

Charged pion candidates are required to be detected in 
at least 12 DCH layers and have a transverse momentum 
greater than lOOMeV/c. Kaons, electrons, and muons are 
separated from pions based on information from the IFR 
and DIRC, energy loss in the SVT and DCH (dE/dx), 
or the ratio of the candidate EMC energy deposition to 
its momentum (E/p). Photon candidates are identified 
with clusters in the EMC with total energy > 30 MeV 
and a shower shape consistent with that expected from 
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a photon. 

The B° J/V^Tr+TT-J^o ^^^^ q- _^ J/^,7r+7r-/i - 
decays are reconstructed in the following way. Elec- 
tron candidates and bremsstrahlung photons satisfying 
2.95 < TO(e+e~(7)) < 3.14GcV/c^ are used to form 
J/tp e"*"e~ candidates. A pair of muon candidates 
within the mass interval 3.06 < m{ii^ ^.~) < 3.14GeV/c^ 
is required for a J/ip — > fJ-'^ n~ candidate. A mass con- 
straint to the nominal J/i/j mass |l3j is imposed in the 
fit of the lepton pairs. We reconstruct K'^ tt+tt^ 
candidates from pairs of oppositely charged tracks form- 
ing a vertex with a probability larger than 0.1%, a 
flight-length significance l/cr{l) > 3 and an invariant mass 
within 15MeV/c2 of the nominal K° mass X(3872) 
candidates are formed by combining J/ip candidates with 
two oppositely charged pion candidates fitted to a com- 
mon vertex. Finally, we form B^{B~) candidates by com- 
bining X(3872) candidates with K^{K^) candidates. To 
suppress continuum background, we select only events 
with a ratio of the second to the zeroth Fox- Wolfram 
moment 0| less than 0.5. 

We use two kinematic variables to identify signal 
events from B decays: the difference between the en- 
ergy of the B candidate and the beam energy, AE = 
E*^ — \fs/2, and the energy-substituted mass mEs = 
^(s/2 + Pj • ^bY/E'I - p|. Here (£'j,Pj) is the four- 
vector (in the laboratory frame) and -^s is the center-of- 
mass (CM) energy of the e+e" system. E'^ is the energy 
of the B candidate in the CM system and ps the mo- 
mentum in the laboratory frame. The signature of signal 
events is Ai? ^ 0, and toes ~ rns where is the mass 
of the B-meson . 

We optimize the signal selection criteria by maximizing 
the ratio n™7(3/2 + ^^) ^ where n™'^ (n™'^) are the 
number of reconstructed Monte Carlo (MC) signal (back- 
ground) events. The optimization was performed by 
varying the selection criteria on AE, toeSi the candidate 
masses of the X(3872) and K'^, and the particle identifi- 
cation (FID) requirements of leptons, pions, and charged 
kaons. The criteria \AE\ < 15MeV, |toes — tob| < 
6MeV/c2 and |m( J/V'tt+tt") - 3872MeV/c2| < GMeV/c^ 
(signal region) were found to be optimal for selecting sig- 
nal events. In case of multiple candidates in an event, we 
select the candidate with the smallest value of |Ai?|. Ap- 
plying our optimized selection criteria, we compute the 
J/V'TT+TT" invariant mass in the range 3.8 — 3.95GeV/c^ 
shown in Figs.^a) and^^b) for the B~ and 5° mode, re- 
spectively. The shaded area shows events in the sideband 
region |toes - 5260| < GMeV/c^. 

We extract the number of signal events with an ex- 
tended unbinned maximum-likelihood fit to the two- 
dimensional distribution y{m-E,Si'mx) where mx is the 
J/'i/'Tr+TT" invariant mass. The probability density func- 
tion (FDF) (normalized to the total number of events) 
is 'P{y) = '^t'^i^tiv) where rit is the number of events 
of category t. We consider three different event cate- 



gories: signal, B decays with the same final-state par- 
ticles as the signal that accumulate near toes ~ ^ns 
(peaking background), and combinatorial background. 
The individual FDFs Vt are assumed to be uncorre- 
lated in toes and mx and can therefore be factorized as 
'Pt{y) = 9t{'mEs)ht{mx), where gt and ht represent the 
Toes and mx probability distributions, respectively. The 
B — > A'(3872)iir signal events are modeled by a Gaussian 
distribution in toes- The resolution function in mx for 
those events is best described by a Cauchy function 
due to the mass constraint of the J/ip candidate. The 
FDF for peaking background events is parameterized by 
a Gaussian distribution in toes and a linear function in 
mx- We model combinatorial background events by an 
ARGUS function |l7l | in toes and a linear function in mx ■ 
The fit performance was validated with MC experiments. 
The mean and width of the toes Gaussian distribution 
for signal and peaking background and the width of the 
mx Cauchy distribution for the B'^ mode were fixed to 
values obtained from MC samples. Other parameters are 
allowed to vary in the fit. 

The fit is performed in the region 5.2 < toes < 
5.3GeV/c2 and 3.80 < mx < 3.95GeV/c2 without 
applying the optimized selection criteria on those two 
variables. The signal region projections of the two- 
dimensional fit are shown in Fig. ^ for the B~ (a,c) and 
B'^ (b,d) modes. We obtain 61.2 ± 15.3 signal events for 
the B~ mode (nj) and 8.3 ±4.5 signal events for the B^ 
mode (n^), respectively. In the following we interpret the 
observed events in the B° mode as the AT (3872). 

The efficiency is determined from MC samples with 
an A:(3872) signal of zero width at 3.872 GeV/c^. The 
decay model consists of the sequential isotropic decays 
B ^ X{3872)K, X{3872) J/ijp°, and p° tt+tt". 
Compared to a three-body decay, this gives a more ac- 
curate description of the observed tt+tt^ invariant mass 
distribution j^. Efficiencies are corrected for the small 
differences in FID and tracking efficiencies that are found 
by comparing data and MC control samples. The final 
efficiencies are (17.4 ± 0.2)% for the B°/K° mode and 
(22.2 ± 0.2)% for the B-fK' mode. 

The branching fraction systematic errors {B^ , B^ 
mode in %) include uncertainties in the number of BB 
events (1.1, 1.1), secondary branching fractions (5.0, 
5.0) 0, efficiency calculation due to limited MC statis- 
tics (0.7, 1.9), MC decay model of the A:(3872) (1.0, 1.6), 
differences between data and MC (1.8, 8.9), FID (5.0, 
5.0), charged particle tracking (6.0, 4.8), and recon- 
struction (-, 1.6). The production ratio of i?" and B~ 
mesons in T(4S') decays is 1.006 ± 0.048 The total 
fractional error obtained by adding the uncertainties in 
quadrature is 9.6% and 12.8% for the B~ and 5° mode, 
respectively. 

Assuming Gaussian systematic errors with a FDF 
Psys{n) ^ exp[-(n - ns)^ /2a'^y^], the negative log- 
likelihood (NLL) function including systematic errors 
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FIG. 1: Signal region projections of m{J/ip -k^ tt ) and jties 
for B- X{3872)K- (a,c) and B° X{3872)K°s (b,d). 
The dashed line represents the combinatorial background 
PDF, the dotted line represents the sum of the combinatorial 
and peaking background PDF, and the solid line the sum of 
all background plus the signal PDF. The shaded area shows 
events in the tties sideband region |mES 
6MeV/c^ 



5260MeV/c^| < 



is Lsys = {[1/L{n)] - [1/lnPsysin)]) ^ where L{n) = 
— hi(C{n) / Crnax) IS the NLL projection of the parameter 
estimate n of the number of signal events and a^ys is the 
systematic error on the number of signal events. The 
significance including systematic errors obtained from 
^y2Lsys(n = 0) is 2.5(7 for the i?° mode and 6.1a for the 
B~ mode. The statistical significance (a^ys = 0) of the 
signal is 2.6a and 7.5ct, respectively. 

Using n|g and rig, the efficiencies, the secondary 
branching fractions and the number of BB events, 
we obtain the branching fractions S° = B{B'^ 
X{3872)K°, X ^ J/i'TT+TT-) = (5.1 ± 2.8 ± 0.7) x IQ-*^ 



and B 



B{B- 



X{3872)K-,X J/iPn+n-) 



(10.1 ±2.5 ±1.0) X IQ-^ For the ratio of branching frac- 
tions, R = B'^/B^, where most of the systematic errors 
cancel, we obtain R = 0.50 ± 0.30 ± 0.05. We calcu- 
late a 90% confidence level (C.L.) likelihood interval 
[ni, Tih] for the number of signal events in the B'^ mode by 
solving the equation 2Lsys{ni^h) = [erf~^(0.95)]^. With 
m = 2.2 and nu = 16.9 the 90% C.L. interval on B" is 
1.34 X 10^6 <B" < 10.3 X 10^6. Using the same strategy, 
the confidence interval on the ratio of branching fractions 
becomes 0.13 < R< 1.10 at 90% C.L. 

We measure the mass of the X(3872) in both modes 
in reference to the precisely measured ipi'^S) mass |l3j . 
We fit the J/i/iTT+TT" invariant mass in the ip{2S) and 



X(3872) region and calculate uix = rnxfit — i^ji!{2S).fit + 
TO^(2s)- The resuh for the mode is (3868.6 ± 1.2 ± 
0.2)MeV/c2 and (3871.3 ± 0.6 ± 0.1)MeV/c2 for the 
B~ mode, where the first error is the statistical un- 
certainty on mx,fit and the second is the uncertainty 
on m^(2S),fit and to^(2S) The mass difference of 

the X(3872) produced in and B~ decays is Am = 
(2.7 ± 1.3 ± 0.2) MeV/c^. The full width at half max- 
imum of the X-mass distribution from the fit on data 
is (6.7 ± 2.7)MeV/c^, which is consistent with the MC- 
determined value of (5.4 ± 0.1) MeV/c^. From this we 
calculate the 90% C.L. upper limit on the natural width 
as r < 4.1MeV/c2. 

Recent observations by BABAR |0| in ISR events pro- 
vide evidence for at least one broad resonance in the 
invariant mass spectrum of J/i^Tr+Tr" at 4.259 GeV/c^ 
that can be characterized by a single resonance with a 
full width of 88MeV/c^. This structure is referred to 
as F(4260). We search in B~ decays for states decay- 
ing into J/ipTT^TT^ above 4GeV/c^ and impose the addi- 



tional selection criterion \m{K' 



-1273MeV/c^| > 



250MeV/c^, which removes backgroimds from 7^1(1270) 
decays. In the resulting mass distribution, Fig. El we 
observe large combinatoric backgrounds and cannot re- 
liably determine the parameters of one or more reso- 
nances. We use a two-dimensional PDF identical to 
the previous model, but fix the central value and width 
of the signal component to the ISR results The 
natural width of 88MeV/c^ has been enlarged by the 
detector resolution, which is found to be the same as 
for the mass region around 3.87GeV/c^. The mx pro- 
jection of the two-dimensional fit is overlaid in Fig. |3 
and yields 128 ± 42 signal events. The statistical sig- 
nificance calculated from \/ —2 ln£o/'C is 3.1ct where C 
and Cq are the maximum likelihood of the fit and the 
null hypothesis fit, respectively. Using a phase-space 
MC simulation of a state at 4.26GeV/c^ decaying into 
J/'i/'TT+TT" and assuming the same systematic uncertain- 
ties and efficiency corrections as for the X(3872), we ob- 
tain the 95% C.L. upper limit on the branching fraction 
By = B[B- y(4260)ii:-,y(4260) ^ J/i/'tt+tt-) < 
2.9 X 10"^ The 90% C.L. hkelihood interval on the 
branching fraction By = (2.0 ± 0.7 ± 0.2) x 10"^ is 
1.2 X 10^5 <By < 2.9 X 10-^ 

In conclusion, our studies of the J/ip-K~^-K~ invariant 
mass below 4 GeV/ yield a signal of 2.5cr and 6. Icr signif- 
icance in the B^ and B^ mode, respectively, with a ratio 
of branching fractions R ~ 0.50±0.30±0.05. We observe 
an excess of events above background in the J/'^tt+tt" in- 
variant mass between 4.2 and 4.4GeV/c^. These events 
are consistent with the broad structure observed in ISR 
events [l]| . 

If one narrow state is observed in the mode B^ 
X{'i%72)K^ , X J/i/)7r+7r^, the diquark-antidiquark 
model Si predicts one amplitude (from Xd or X^) to be 
dominant in the charged mode and the other amplitude 
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FIG. 2: Projection of tlie two-dimensional fit for the mass 
region above 4GeV/c^ for events within the tties signal region 
in B~ —> J/iIj TT^TT^ K~ . The lines represent the same event 
types as in Fig. 



to be dominant in the neutral mode. In this case, the 
model predicts the relative rates to be equal (i? = 1) and 
the mass difference to be (7 ± 2)MeV/c^. The ratio of 
branching fractions is consistent with our measurement, 
0.13 < R < 1.10 at 90% C.L., and the observed mass 
difference of (2.7 ± 1.3 ± 0.2) MeV/c^ is both consistent 
with zero and the model prediction within two standard 
deviations. In the S-wave molecule model 0, the neu- 
tral mode branching fraction is predicted to be at least 
10 times smaller {R < 0.1) than the charged mode. How- 
ever, we obtain a ratio of neutral to charged branching 
fractions which is slightly more consistent with isospin- 
conserving decays. 
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